Semiclassical analysis of edge state energies in the integer quantum Hall effect 
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Analysis ol edge-state energies in the integer quantum Hall effect is carried out within the semi- 
classical approximation. When the system is wide so that each edge can be considered separately, 
this problem is equivalent to that of a one dimensional harmonic oscillator centered at x = x c and 
an infinite wall at x = 0, and appears in numerous physical contexts. The eigenvalues E n (x c ) for 
a given quantum number n are solutions of the equation S(E,x c ) — n[n + j(E,x c )] where S is 
the WKB action and < 7 < 1 encodes all the information on the connection procedure at the 
turning points. A careful implication of the WKB connection formulae results in an excellent ap- 
proximation to the exact energy eigenvalues. The dependence of "/[E n (x c ), x c ] = jn(x c ) on x c is 
analyzed between its two extreme values | as x c —* —00 far inside the sample and § as x c — > 00 far 
outside the sample. The edge-state energies E„(x c ) obey an almost exact scaling law of the form 
E n (x c ) = 4[n + 7n(zc)]/( ^4^ +3 ) an d the scaling function f(y) is explicitly elucidated. 

PACS numbers: 73.43.Cd, 03.65.Sq 



I. INTRODUCTION AND STATEMENT OF 
THE PROBLEM 

The concept of edge states is central for elucidating 
the physics of quantum Hall and spin quantum Hall 
systems 1 -^^. It is illustrated by considering an elec- 
tron (mass m and charge — e) restricted in two dimen- 
sions to a stripe —00 < y < 00, — L x < x < 0, and 
acted upon by a perpendicular magnetic field B = Bz = 
VxA. At strong enough magnetic field it is reason- 
able to assume L x ^> I, (where I = \J~^ is the mag- 
netic length), and consider a system with a single edge 
at x = 0, where the electron is confined on the half plane 
— 00 < y < 00, — 00 < x < 0. Within the Landau 
gauge A = (0, — Bx), the Schrodinger equation and the 
corresponding boundary conditions read 



2^ 



(p y + -Bx) 2 ]^{x,y) 



¥(0,y) = ¥(-oo,y) = 0. 



E9(x,y), (la) 
(lb) 



Translation invariance along y enables the replacement 
p y — ► k and turns equation {Taj into that of a one di- 
mensional oscillator centered at X c = kl 2 with hard wall 
boundary condition at x = 0. Using I as unit of length, 
and huj c /2 = ti 2 /2ml 2 as unit of energy, the correspond- 
ing eigenvalue problem then reads (—00 < x < 0) : 



d 2 , 
'dx^ + {Xc 



x) 2 ]^{x) 
^i(O) = V(-oo) = 0, 



(2a) 
(2b) 



where all quantities are dimensionless and — 00 < x c = 
X c /l = kl < 00 is the (dimensionless) center of the oscil- 
lator, which is considered as a continuous parameter. In 



the absence of edge, the eigenvalues are simply the Lan- 
dau energies E n = 2n + 1, independent of x c . The pres- 
ence of edge changes this infinite degeneracy and turns 
the energies to be dependent on x c . When \x c \ is small, 
the oscillator center is close to the edge within a mag- 
netic length so that the corresponding eigenfunctions are 
localized near the edge of the system at x = and hence 
they are referred to as edge states. 

The challenge of calculating the eigenstates ip n {x;x c ) 
and the eigenvalues E n (x c ) is referred to as (partially) 
restricted harmonic oscillator problem, and dates back 
several decades ago^. Since confined quantum mechan- 
ical systems are ubiquitous, it arises in numerous con- 
texts in physics and chemistry, much beyond the specific 
edge-state scenario of the quantum Hall effect mentioned 
above. Formally, the eigenvalue problem (|2al2b|) pos- 
sesses an exact solution in terms of Kummer functions 6 , 
but this solution is hardly useful, as these functions are 
rather complicated. 

The present work focuses on the eigenvalues E n {x c ) 
of the Schrodinger equation (|2a|). analyzing their depen- 
dence on x c and the energy quantum number n. It 
has two main objectives. The first one is to examine 
the semiclassical approach for estimating the eigenvalues 
E n (x c ). The idea of employing the WKB approxima- 
tion for the study of confined systems has a long history. 
Probably, the discussion most relevant for the present 
study is a four decade old work by Vawter— , but see 
also references ^ 10 ' 11 ! 12 ' 13 ' 14 . Here the WKB method is 
augmented and adapted to the special problem at hand. 
The main subtleties resulting from the fact that there are 
three regions in [x c , E) space (where the WKB connec- 
tion formulae assume different forms) are clarified, and 
the crossover between these regions is treated with special 
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care. The energy eigenvalues are found from an expres- 
sion S{E 1 x c ) = 7r[n + j(E, x c )] relating the WKB ac- 
tion S(E, x c ) at energy E with a function j(E, x c ) which 
encodes all the physics behind the connection formulae. 
The main achievement of this part is an elucidation of 
"f{E, x c ) for the entire range E > 0, — oo < x c < oo and 
a subsequent evaluation of the eigenvalues E n (x c ) which 
appear to be in excellent agreement with the exact ones. 
As a byproduct, the source of disagreement (as noticed 
in Ref»I) between the WKB solutions and the exact en- 
ergies occurring when the right turning point is close to 
the edge, is remedied. 

Our second goal is to examine whether a scaling re- 
lation of the form E n /L" = f(x c /Ln^ 2 ), exists within a 
wide domain of energies E n and oscillator center parame- 
ter x c . Here L n is a linear function of the energy quantum 
number n = 0,1,2,... which serves as a length scale in 
energy space, while a is a scaling exponent. Our main 
result in this context is that there is indeed an approx- 
imate scaling law (albeit with small correction), which 
takes the form 



E n (x c ) 
4n + 3 



/( 



3 - 4j n (x c ) 
An + 3 J 



(3) 



The function f(y) is a solution of a simple implicit equa- 
tion , (the scaling variable y = x c f\J An + 3 should not 
be confused with the y coordinate of the particle) . The 
scaling relation ^ of edge state energies is supported by 
numerical results based on exact diagonalization. 

The rest of the paper is structured as follows: In sec- 
tion |TT] the principles of the WKB method as employed 
here are recalled in terms of the WKB action and the 
function j(E^x c ). It becomes evident that the magni- 
tude and sign of well as those of the right turning 
point crucially determine the way of constructing the con- 
nection formula. The relevant algorithm is worked out 
and and the WKB approximated energies E n (x c ) are ob- 
tained and compared with the exact ones in section HTT1 
Finally, in section HVl the scaling hypothesis is tested and 
a scaling relation for E n (x c ) is suggested. 



II. SEMICLASSICAL QUANTIZATION 

This section briefly introduces and discusses the ba- 
sic ingredients of semiclassical action and the function 
j(E,x c ) which are necessary for the subsequent imple- 
mentation of the WKB analysis. For notational conve- 



nience the harmonic term (x 



) 2 with the hard- wall 



condition at x = are combined into a single potential, 



v(x) = ( ' T Xc)2 ' 



x < 
x > 



(4) 



In the WKB method applied for a bound-state problems, 
the energy E is higher than the potential energy v(x) 
between the two classical turning points x± < X2 beyond 



which the solution is classically forbidden. A glance at 
figure Q] indicates that the turning points are, 



Xi 



x 2 = min(a; c + Ve, 0), (5) 





Xc 



FIG. 1: Harmonic oscillator potential (x — x c ) 2 and a confining 
wall (edge) at x = 0. Several regions in (E, x c ) parameter space 
have to be considered. For x c < 0, the low energy region (/) is 
defined by E < x 2 , and the high energy region (Ila) is determined 
by E > x 2 . For x c > 0, there is only a high energy region lib for 
which E > x 2 c . The crossover region is marked by the horizontal 
dashed line in the left panel. 

The following definitions and notations will be used 
below in discussing some general features, 



k(x) = ^JE - v(x) > 0, (x 1 < x < x 2 ), (6a) 
S(x, x') = ( k(t)dt, (x t < x < x < x 2 < 0X,6b) 

J X 

S(E,x c ) = S(xi,x 2 ). (6c) 

The dependence of the WKB action S(E, x c ) on energy 
enters through that of the turning points as in equation 
([5]). This action should satisfy a certain relation result- 
ing from the WKB connection formulae at the turning 
points X\ < x 2 . In most cases, the structure of this 
relation involves the action S{E 1 x c ) as an argument of 
a trigonometric function, and its decoding can then be 
casted in the general form, 



S(E,x c )=7r[n + i(E,x c )]. 



(J) 



It should be emphasized that equation Q is not an iden- 
tity, but, rather, an implicit equation for finding the 
eigenvalues E n {x c ). The function < ^{E^Xc) < 1 
should be calculated independently and encodes the de- 
tails of the connection formulae. Arriving at this equa- 
tion (equivalently knowing ^(E, x c ) explicitly) is the cen- 
tral part of the WKB method (the computation of the 
action S(E,x c ) is straightforward in most cases). Once 
this task is achieved, the solutions E n (x c ) of this implicit 
equation are the required eigenvalues within the WKB 
approximation. Inserting the solutions into equation ([7]) 
implies the identity 



S(E n (x c ),x c ) = 7r[n + ~f(E n (x c ),x c )} 



(8) 



which defines ~/(E n (x c ), x c ) = r y n (x c ) self consistently as 
function of the discrete energy quantum number n and 
the continuous variable x c . 

The easier part of the WKB procedure is calculation 
of the action itself, which in the present problem of con- 
fined harmonic oscillator becomes elementary. It depends 
on the appropriate region in (E, x c ) parameter space, as 
depicted on figure [TJ In the low energy region, / the 
positions of the turning points are x\ — x c — \f~E~ and 
X2 = x c + \f~E~ and the action S is given by 



S = 



-Ve 



1 



r — x c ) 2 dx — 2 n E- (9) 



In both high energy region Ila, lib, the turning points 
are located at x\ = x c — \f~E~ and X2 = 0. The action is 
now dependent on x c and it is given by 



n E/2 




FIG. 2: The action S(E, x c ) = ^s(t) as a function of t = x c /\/E 
(with E fixed). It varies between n E/2 for t = -1 (ie c < -s/E) 
and for t = 1 (x c > 0). For t = x c = 0, it is equal to 7r£/4. 







Ve 
ttE E 



2 dx 



arcsm 



E 



■VE- 



(10) 



The action has a dimension of energy and can be written 
as 



S(E, x c 



ttE 



s(x c /VE), 



where the dimensionless function 
1 1 



s(t) = — arcsint \J\ — t 2 

2 7T 7T 



(11) 



(12) 



has the following expansions 



t ~ ±1 
t ~ 



-(*) = ^ ± ^(1 - l*l) 3/2 (13a) 
s(t) = \- 2 4- (13 b ) 



The dependence of the action S(E,x c ) — ^s{t) on t = 



x c /vE is displayed in figure 

The less trivial task which will occupy us in the next 
subsection is to calculate the function j(E,x c ). When 
the potential is smooth and the wall is far from the right 
turning point X2, 7 is determined by the standard WKB 
procedure such that the solutions for x < x\ and x > xi 
(the classically forbidden regions) contain only decaying 
exponents. Each turning point Xi then contributes a term 
^ where /3, is an integer (also referred to as the Maslov 
index for the pertinent turning point), and 7 = + 
ftz). This is the situation appropriate in the lower part of 
region / in figure [1] where (3\ = fa = 1, implying "f = \- 
The definition can be extended to the situation when a 
turning point occurs at the wall. Although there is no 
solution to match on the right of the turning point, the 
condition ip(0) = replaces the connection formula and 



implies P W al\ = ^" This is the situation relevant for the 
higher part of region I la and for region lib in figure [TJ 
implying 7=1- As will be stressed below, the situation 
is different in the crossover region (around the horizontal 
dashed line in figure [T|). Since j(E,x c ) is expected to 
vary smoothly between i and f , it cannot be written 



in terms of an integer fa. 
parametrization 7 = \(J5\ 



If one insists on the same 
fa}, this implies that fa is 



a non-integer Maslov inde»i^. An attempt to use 7=5 
for E < x 2 and 7 = 

to an artificial discontinuity at E = x 2 since 

S(E = x 2 c + 0,x c ) = ^=n(n+~ 



I for E > x 2 is too naive and leads 
,2 



— > E n (x c ) = 2n+ -, 

S(E = x 2 c - 0, x c ) = 
— > E n (x c ) = 2n+l 



irE , 3, 
2 V 4' 



(14a) 
(14b) 



This artificial discontinuity is displayed in figure [3] below. 
As will be explicitly demonstrated in the next section, the 
correct picture is indeed different: for fixed n, 7„(a; c ) is 
a monotonic and smooth function of x c which tends very 
quickly to \ at negative x c and to | at x c > 0. 



III. WKB EVALUATION OF EIGENVALUES 

A proper analysis of the crossover regions within the 
WKB formalism should then remedy the discontinuity 
problem by a proper treatment of the medium energy re- 
gion E ~ x 2 . This is carried below by approaching the 
cross over region from below (E <C x 2 — > E w x 2 ) and 
from above (E ^> x 2 — * E » x 2 ). The case x c ~ where 
the oscillator center is very close to the wall requires some 
special treatment at low energy where the linear approx- 
imation for the potential near the turning point is not 
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FIG. 3: Energy eigenvalues E n (x c ) obtained within the semiclas- 
sical approximation J7), where S[E, x c ) is given in equations I|9ll0j l 
while 7 = \ for E < x\ and 7 = | for E > x\. The (artificial) 
discontinuity occurring along the dashed line E = x\ results from 
an improper treatment of the connection formula in the crossover 
region. 



valid. After this task is completed, the eigenvalues cal- 
culated by the WKB approach are compared with the 
exact ones, and the agreement is virtually perfect. 



A. Approaching E ~ x c from below 



When E < x\ (region / in figure ([T])), the turning 
points and the WKB action are 

£1,2 = x c ± vE, (15a) 
S(E) = 1 --kE. (15b) 

For E <C x\ the eigenvalues are expected to approach 
the energies 2n + 1 of the unrestricted harmonic oscil- 
lator. The question is how they are modified when E 
approaches x\ from below. For energies close to x c the 
turning point X2 is rather close to the wall at x = 0, and 
this must be taken into account. Practically, it means 
that the wave function in the classically forbidden region 
X2 < x < must include also an exponentially increasing 
contribution (beside the exponentially decreasing one), 
since vanishing of the wave function at the wall x — 
can be achieved only by a proper combination of the two 
waves. Therefore, a modification of the standard WKB 
connection procedure is required near X2 in order to take 
into account the effect of the wall behind the turning 
point. On the other hand, the connecting procedure at 
the left turning point is standard. Recall that to the 
left of the point x\ there is a single wave, that decays as 
x — > — 00. This means that within the linear approxima- 
tion for the potential near the turning points, only the 
first Airy function Ai is used near X\. The connecting 



procedure around x\ implies that the WKB approxima- 
tion for the wave function for x\ < x < X2 < is 



ip(x) 



C 



C 



: sin[<S(a^i, x) + 



VKx] 



sm[S — S(x, X2) 



7T. 



(16) 



where C is a constant. The modification near the second 
turning point X2 which is close to the wall assumes that 
in the small region X2 < x < the potential is approxi- 
mated by a linear expansion near X2 ■ 

(17) 



v(x) w E + 2VE(x - x 2 )- 
In terms of the variable 

z = z(x) = (4E)i(x ~ x 2 ), 



(18) 



and within the approximation (|17[) . the Schrodinger 
equation at the vicinity of X2 reads 



d 2 iP(z) 
dz 2 



zip(z) = 0. 



(19) 



Unlike the analysis near X\, here the combination of both 
Airy functions is required in order to construct the solu- 
tion around X2, that is, 



ip(z) = ocM(z) + /3Bi(z). 



(20) 



When this combination is examined at z < (that is, 
x < X2) and the asymptotic form of the Airy functions is 
used, it is found that the wave function at x < a; 2 is, 

ip{x) = asm[S(x,x 2 ) + -] + (3 cos[S(x,x 2 ) + — 



Since this is the same function as that found in (fTBj) we 
can compare the two and obtain a constraint on the co- 
efficients C,a,fl. To this end, denote ( = S(x,X2) + f , 
and open the sine and cosine functions, to get 

sin C(C* sin S + a) + cos ((C cos S + (3) = 0. (22) 

This holds for every £, in the appropriate interval, hence 
each term in the bracket should vanish separately. This 
yields a couple of equations for the three coefficients 
C,a,f3. The wave function (|20[) should vanish at the 
wall x — 0, that is (zq = z(x — 0)) : 

ip{z ) - 0, (23a) 

z = -(4E)*x 2 = -(4E)^(x c + VE)>0, (23b) 

implying that 



R A1 (^0) 

Bi(z ) 



(24) 



Employing the statement after equation ([2T]) this leaves 
us with two homogeneous equations for C and a : 



CsinS"- 
CcosS 



a = 0, 
Ai(z ) 
' Bi(z )' 



(25a) 
(25b) 
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and the energies E n are obtained by requiring the van- 
ishing of the determinant 



cos S 



Ai(z ) 
Bi(z ) 



sinS 1 = 0. 



(26) 



The implicit dependence on the energy E enters through 
the dependence of S and Zq on energy, equations (|15b[) 
and (123bj) respectively. Equation (|26|) is recast in the 
WKB format as 



ttE 



S = '—r- = n[n + j(E, x c )] 
where < 7 < 1 is explicitly given by 



j{E,x c ) = 1 arctan — — -. 

7T Al(z ) 



(27) 



(28) 



with z defined in (|23b[) . The function j(E, x c ) is plotted 
as function of the energy E for several values of x c in fig- 
ure [H The eigenvalues E n (x c ) are the intersection of the 



This value 7 = | is intermediate between the "high" and 
"low" energy values | and \. For n = this gives the 
energy | instead of | and 1. The exact value is indeed |. 
The disagreement between the "high energy" prediction 
I and the exact value | has been noticed for n = in 
table II of Ref,» where it is dubbed as a WKB error. 
Our results show, however, that the WKB formalism as 
applied here is virtually exact. 



B. Approaching E ~ x\ from above 

Consider now region (II) in figure [TJ where E > x 2 . In 
order to treat this "high energy" region the linear expan- 
sion of v(x), eq. (fTTf is now replaced by 



v(x) ~ x 2 — 2x c x, 



(32) 



so that the wave function satisfies the Schrodinger equa- 
tion and the hard wall boundary conditions at a; = 0, 
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FIG. 4: The function ^(E,x c ) defined in equation l|28j l is plotted 
versus E for different values of x c = —4, —3, — 2, — 1, 0, 1, 2. 

curves j(E, x c ) (as function of E with Xq cLS S- parameter) 
with the straight lines —S(E, x c ) — n. In the two limiting 
cases and E w x 2 these solutions can be easily 

evaluated. When E <C x 2 , the value of zq is large. Then 

1 3 

(29) 



Ai ( 2 o) _ ~ 
2 



Bi(z ) 



3 "0 



0. 



and equation (f28[) leads to the familiar result of the har- 
monic oscillator problem with 7=1/2 and E n = 2n+ 1. 
On the other hand, when the energy is very close to (but 
still lower than) x 2 then zq is very small, and hence, 



Ai(z ) _ Ai(0) 
Bi(z ) ~ Bi(0) 

so that 7 = 2/3 and 

E n = 2n + 



1 

V3' 



(30) 



(31) 



ip" - (x 2 c - E - 2xx e )ip = 0. 
tp(0) = 0. 



(33a) 
(33b) 



whose solution has the form 

x r ' " _£/ — — ^xxn 



tp(x) — aAi 



|2x,l 2 /3 



/3Bi 



|2x c | 2 /3 



(34) 



The hard wall condition at 2; = implies 



a 



Ai 




|2£C C | 2 / 3 


Bi 


" X 2 r -E ' 


\2x c \ 2 /3 



_ Ai(-W) 
= Bi(-W) : 



where 



(35) 
(36) 

(37) 



At large distance x < 0, the asymptotic form of the 
Airy functions for negative arguments can be used to 
yield 



ip(x) oc a sin 
+ /3cos 



r x 2 — E — 2xx c 
x r , _£/ ^iXXf 



+ 



|2:r c | 2 / 3 ' 4_ 
and the action for the linearized potential v(x) is 



(38) 



«S(z,0) 



J X 



.i r* dec 

(E - x 2 c f/ 2 (E-x 2 c + 2xx c f' 2 



SXr 



3x c 



.(39) 
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Therefore the wavefunction has the form 



i/j(x) oc a sin — cS(x, 0) + 6 + — 



+ /3cos 



-eS(x,0) + 5 + 



with 



e = sign(x c ), 



5 = 



3\x c \ 



3 



(40) 

(41a) 
(41b) 



On the other hand, following the procedure leading to 
equation (|16[) the wave function resulting from the con- 
necting formula at x — x\ is (x% < x < 0) 

C TT 

ip(x) = — t== sm.[S — S(x, 0) 



-1=00^0)^-5]. 



(42) 



It is now allowed to compare the two expressions (|40[) 
and (|42[) for the same wave function. These two expres- 
sions should be equal for any value of the variable S(x, 0). 
Equating the coefficients of sm<S(a:,0) and cosS(x,0) to 
(each one separately), and employing the relation (|3"5)) 
then result in the following two equations for C and a : 

- aecos(<5 + ^) + /3esin(<5 + y) = -Ccos(S + 



asm(5 + — ) + (3cos(6 + — ) 



Csin(S'+ -X43a) 



leading to 



with 



tan(S + -) = e cot(<5 - 5 + -) 



a _ Bi(-W) 



tan(5 = — = 

P 



Ai(-W) ' 



(44) 



(45) 



Finally, the energy levels E(x c ) are given by the action 
M and 



,r, *\ ( m~W\ 2 Tir3/2 TT 

tan(S + -) = -e cot arctan A . tt/ + -W 3/2 + - 
4 y V Ai[-W] 3 4, 

(46) 



or, equivalently, by the central expression, 
S(E, x c ) = Tr[n + j(E, x c )] 



where < j(E,x c ) < 1 is given explicitly by (compare 
with equation (|28|) applicable for the case E < x 2 ), 

l(E,x c ) = \ 

e ( Bi 2 TT 3 tt\ 

arctan cot arctan- — \- - W 2 H — (48) 

7r V Ai 3 4 / 



This function, calculated here for x c > —\JE is shown on 
figure ([4}. It has simple values in the following special 
cases : 

x c = -VE, W = 0, j(E,x c ) = 2/3, E n = 2n+^, 

(49a) 

x c = 0, W = oo, i(E, x c ) = 3/4, £„ = 4n + 3. 

(49b) 



C. Perturbation expansion for small \x c \ 

Near x c — and at low energy the linear approxima- 
tion is not justified because at x c — the linear term 
vanishes. It is then necessary to modify the semiclassical 
formalism for evaluating E n (x c ) and 7 n (^c) near x c = 0. 
Here a simple perturbative expansion for small x c is used. 
The unperturbed hamiltonian is, 



H Q =p 2 + x 2 



(50) 



The solutions tp n (x) with ip(0) = are the antisymmet- 
ric eigenfunctions of the free harmonic oscillator (with 
energy E n = An + 3), 



e-* 2 / 2 H 2n+1 (x) 
x / 7 r2 2 "(2n+ 1)!' 



(51) 



where H n are the Hermite polinomyals. The perturba- 
tion term for small and the first order cor- 
rection then yields the perturbed energies, 

,o 

E n (x c ) w 4n + 3 — 2x c I xip^(x)dx 

J — oo 

= ( 4 » +3 ) + ^n( i+ 5) (52) 



The value of 7 c (x c ) near x c = can therefore be obtained 
from the quantization condition eq. |8[). Expanding the 
action S(E,x c ) near a; c = to first order in x c reads 
S(E,x c ) — nE/4 — \fEx c , which implies, 



7n(x c ) = ^ + /? n a; c . 



(53) 



(47) The coefficients /3 n are given by, 



0o = 0.01286 

#l = 0.00416 

f3 2 = 0.00214 

/3 3 = 0.00135 



V4n + 3 



(54) 



7 



In the limit n — > oo, the product tends to y 3 , so that 

Poo = 0. Equation (|53p shows that at low energy and for 
x c very close to (when the linear approximation to the 
potential fails), ^ n {x c ) shoots up slightly above §. Since 
f3 n is small, the deviation (3 n x c is virtually negligible. 



A. Scaling of y„(x c ) 

Before discussing the scaling of the energies themselves 
it is useful to check a possible scaling behavior of 7 n (x c ). 
To obtain this function explicitly the eigenvalues E n (x c ) 



D. Comparing WKB results with the exact ones 

The central result of the foregoing discussion can now 
be presented. Once the function j(E, x c ) is known, the 
spectrum E n (x c ) is entirely determined by the implicit 
equation ([7] ) . This is calculated and displayed as solid 
lines in figure ([5|). The exact eigenvalues are marked by 
dots on the same figure. The fit is indeed perfect. 



TP 
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FIG. 5: The solid lines display the energy levels E n (x c ) calculated 
in the WKB approximation as function of the continuous parameter 
x c . The exact eigenvalues (obtained by numerical diagonalization) 
are represented by the full circles. 
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FIG. 6: The function "fn(xc) as function of x c for n = 0,1,2,3. 
Its calculation is explained at the beginning of this subsection. 

are calculated as explained in the previous section, and 
then substituted instead of E in the function ^(E, x c ) 
defined in equations (|28I48[) ) leading to the desired func- 
tion j n (x c ) = 7[EVj(x c ), x c ]. It is displayed as function 
of x c for several values of n in figure © . Looking at fig- 
ure [6] one is tempted to search for a " scaling relation" in 
the sense that "/ n (zc) depends on a certain combination 
X{n, x c ) of x c and n. Such combination can then be used 
as a scaling variable where all the curves collapse on a 
single one. Based on our previous analysis it is estab- 
lished that J n (xc) an d the spectrum E n (x c ) is very well 
fitted (within 2 %) by the function 

ln{x c ) = 1 (55) 

x Bi[-2 1 /3( 2 „+ l)^(x c + J2n+ §)] 

— — arctan — 

n Ai[-2V3(2n + |)V 6 (a; c + J2n+±)] 



IV. INSPECTING SCALING BEHAVIOR for x c < ~ V 2n + f , and by the function 



After substantiating the efficiency of the WKB method 
for the confined harmonic oscillator problem (summa- 
rized in figure O, and elucidating the functional form 
of j(E,x c ) (see equations (|28l48p ) let us return to the 
question of scaling posed in the Introduction. It should 
be stressed that the notion of scaling here simply means 
that some functions of x c and n can be represented as 
functions of a certain combination of these two variables. 
It has nothing to do with thermodynamics, of course, al- 
though the variable n is, in some sense, analogous to the 
notion of length scale. 



, x , 1 f Bi[-WJ 2 o /2 \ 

ln{x c ) = 1 arctan tan arctan + - W*' 

7T \ Ai[— W„J 6 ) 



with W n (x c ) = 2n+ - -^|2x c | 2 / 3 , 



(56) 



for — \/2n + | < x c < and by 7 n (x c ) = 3/4 for x c < 0. 



Strictly speaking then, there is no single combina- 
tion of n and x c that enters both expressions. Practi- 
cally, however, when the variable X = (2n + |) 1 / 6 (x c + 



2n + |) is used, the collapse of all curves is very good 

as can be seen in figure [7] where 7 n (x c ) = j{X), is plotted 
against X. 




FIG. 7: The functions "f n (x c ) are plotted as a function of the 
scaling variable X = (2n + |) 1 / 6 (a' c + ^2n + |). The dashed line 
is the function 'y(X) defined by equation l|57|l . 



Moreover, ^(X) is very accurately fitted by the func- 
tion 



7 (X) 



1 



1 



4 2expAX + I 



(57) 



where A ~ 3.5, represented by the dashed line in figure 
[7J The energies E n (x c ) are then calculable (within 2 
% accuracy) through the solution of the standard WKB 
implicit equation, 

S{E n ,x c )=ix[n + 1 (X)] (58) 

with 'y(X) as given by the approximate formula (|57l) . 

B. Scaling of edge state energies E n (x c ) 

Turning now to the possible scaling of edge state ener- 
gies, consider the identity (JSj) and suppose, for the mo- 
ment that "fn(xc) is a constant. After dividing both sides 
by 27r(n + 7) and defining the scaling variables 



y = x c /2^n + 7, 
u = E/4(n + j), 



(59a) 
(59b) 



the identity ([8|) takes the form 



1 — u u U V 1 7 

— 1 — arcsin^=H — \/u — y 2 = 0. (60) 

2 7T sJU 7T 

Equation (|6T)|) then implicitly defines the functional rela- 
tion 



En 



4(n + 7) 



= m = /( 



), (61) 



f(y) 
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FIG. 8: Universal function f(y) and y 1 (dashed line). 



between u and y. The function f(y) is shown in figure [5] 
below. From the expansions (|13b[) of the function s(t), it 
is easy to deduce that the function f(y), shown on figure 
([!]), has the following limits 



1 

2' 



V < -1/V2 f(y) 



0>y> -1/V2 f(y) = \ + + y) 3/2 , 



4y 



y~0 /(j,) = l + ^, 

7T 



2 VV2 



2/3 



., 4 / 3 



(62) 



However, the relation (|6H)) is of little use, because 7 is 
not known in advance and its evaluation is part of the 
problem. A more practical approach would then be to 
fix an appropriate constant value 70 and take account 
of the small variation of | < 7 < | by some correction. 
From the edge-state point of view, the region x c > — 1 for 
which 7 I = 70 is the most relevant one. Therefore, 
the "length scale" is fixed as L n — 4(n + 7o) = 4n + 3 and 
expand equation (|6ip to first order around 70. To save 
notations u and y are now defined with 7 = 70 = 3/4. 
The result is 



E„ 



4n 



f(y) 1 - 



3-47 

4n + 3 



1 - 



yf(y) 
2/(2/) 



(63) 



Strictly speaking, the problem mentioned in connection 
with equation (|60|) still remains, but this time it is casted 
in a form showing that exact scaling is expected to work 
for large n as L n = 4n + 3 — > 00. The second term in 
the square brackets can be neglected because for large 
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2/7 = 3/4 and the correction disappears, whereas for 
small y, f(y) — * 1/2 so its derivative vanishes. It is then 
reasonable to suggest the following relation 



E n (x c ) 



= m, 



4[n + 7„(x c )] 
where, employing the expression 157)) for ~/ n ( x c), 

1 



4{n + ln (x c )} = 4n + 3 



2 e AX Tl (x a ) + i> 



(64) 



(65) 



with X n {x c ) = (2n+^) 1 / 6 



(2n- 



with A =~ 3.5. To test the scaling relation, the exact 
eigenvalues E n {x c ) for n — 0, 1, . . . 13 were calculated 
via numerical diagonalization resulting in 14 curves each 
contains 100 points x c (we have already shown that these 
energies can also be calculated within the WKB formal- 
ism developed here). The 1400 numbers appearing on 



-3 



-2 
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FIG. 9: Test of the scaling relation H64I I. The numbers E n (x c ) 
are obtained numerically by exact diagonalization and the func- 
tion 7(a; c ) is approximated employing equation {65) . The numbers 
appearing on the LHS of equation i1(i4[) are then displayed as func- 
tion of the scaling variable y = ,f c „ , thus generating a smooth 
curve representing the scaling function f(y) given by eqs. (160161ft . 
It results from the collapse of 14 curves E n (x c ) with n = 0, 1, . . . 13, 
each contains 100 points x c . 



the LHS of equation ([64]) are thus computed and dis- 
played against the scaling variable y in figure [5] below. 
The smooth curve gives the scaling function f(y) and for 
y > — -i= it coincides with the curve of figure [3 The 
collapse of 14 curves on a single one is remarkable. The 
upshot is that equation (|64p with f(y) analyzed in equa- 
tions (|6"2")) and ^{x c ) parametrized as in equation ([57)) is 
an excellent approximation for E n (x c ). 
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